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Spontaneous parametric down conversion has been shown to be a reliable source of entangled 
photons. Amongst the wide range of properties that have been shown to be entangled, it is the 
orbital angular momentum that is the focus of our study. We investigate, in particular, the bi- 
photon state generated using a Gaussian pump beam. We derive an expression for the simultaneous 
correlations in the orbital angular momentum (£) and radial momentum (p) of the down-converted 
Laguerre- Gaussian beams. Our result allows us, for example, to calculate the spiral bandwidth with 
no restriction on the geometry of the beams: £, p and the beam widths are all free parameters. 
Moreover, we show that, with the usual paraxial and collinear approximations, a fully analytic 
expression for the correlations can be derived. 



PACS numbers: 42.50.Tx, 42.50.Dv 

INTRODUCTION 

Entanglement is one of the most remarkable phenom- 
ena introduced by quantum physics. It has applications 
in quantum imaging []]-[4] and quantum cryptographic 
systems [5[ and is of great interest both within the con- 
text of quantum information theory [6[ and in quan- 
tum computing 0, H| . More recently, entanglement has 
been demonstrated between spatial modes carrying or- 
bital angular momentum (OAM) and this has been used 
for quantum information protocols [9l4l2| and to demon- 
strate violation of Bell's inequality [l3|, [l4| and in an 
angular demonstration of the EPR paradox [IBj]. En- 
tanglement has also been shown in the spatiotemporal 
structure of the light, both for type-I [l6[ and for type-II 
[Tt( | parametric down-conversion. 

Spatial modes can also carry radial momentum. Both 
angular (£) and radial (p) momenta can be controlled and 
converted between one another using diffractive optical 
elements tipically in the form of spatial light modula- 
tors (SLMs) displaying holograms [18|, fl9j . Recent ap- 
plications of controlling the entanglement in the radial 
modes involve topological knots of phase singularities. 
These can be produced by engineering specific superpo- 
sitions of radial and orbital angular momentum modes 
[l8| . Topological knots have been shown to be robust 
against perturbations [20], a feature that might be help- 
ful when propagating entangled states through the atmo- 
sphere. Another application of the entanglement in the 
radial modes is a larger number of qubits that could be 
used in a quantum information task. 

In spontaneous parametric down-conversion (SPDC), a 
pair of lower- frequency photons are created when a pump 
field interacts with a nonlinear crystal [9|, [2l| . The spatial 
structure of the down- converted photon states can be ex- 
pressed as a mode decomposition of their joint wave func- 
tion in an appropriate basis. The degree of entanglement 
is directly related to the number of modes participating 
in the state, often referred to as the spiral bandwidth 



[22I , [23| . The strength of the entanglement between pairs 
of photons can be verified in the laboratory by converting 
with SLMs specific £ modes to the fundamental mode of 
the single mode detection fibers and measuring the coin- 
cidence rate [20|. 

For applications in quantum information, such as mul- 
tidimensional quantum imaging 0, [HI, or f° r applica- 
tions involving SLMs, which affect both the £ and the 
p quantum numbers jl9|, [25| , it is important to be able 
to calculate the exact form of the down- converted bipho- 
tons, with no restrictions on any of the denning parame- 
ters: the OAM indices 1^ £ S: the radial indices p^, p s or 
the beam widths w p , Wi, w$ . 

In this paper we consider a typical SPDC set-up as 
shown in FigJTJ We first analyze the general problem, 
making use of the usual assumptions for the phase match- 
ing (i.e. we assume that the transverse phase matching 
is satisfied and that there is no extra longitudinal phase 
mismatch, other than the difference of the k vectors of 
the down- converted photons) and derive an integral so- 
lution for the mode amplitudes. This is simple enough to 
evaluate numerically and we show that it reduces to an 
analytical solution in the (experimentally relevant) limit 
of a short crystal. Special cases of our results are in ex- 
cellent agreement with previous calculations (23|, when 
the relevant restrictions are applied. We show, moreover, 
that the same analytical result is obtained, without the 
need for a thin crystal approximation, when we consider 
collinear geometry. 

We use our results to investigate the influence of the 
relative beam widths and the effect of the radial contri- 
butions on the resulting spiral bandwidth. We show that 
the various mode amplitudes have a strong dependence 
on the ratios of the signal and idler widths to the pump 
width. In a future paper we will discuss the effect of 
the crystal length, phase matching, and the differences 
between collinear and non-collinear geometries. 
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GENERAL AMPLITUDES OF THE SPDC STATE 
Geometry and notation 

A typical SPDC setup consists of a continuous-wave 
Gaussian pump beam propagating in the z direction, as 
shown in FigJTJ incident on a short (typically l-3mm) 
nonlinear crystal of length L. This produces two highly- 
correlated, lower-frequency photons, commonly termed 
signal and idler. Energy is conserved in this process so 
that uj p = u s + cj$, where the subscripts p,s,i refer to the 
pump, signal and idler, respectively The photons are 
emitted at angles to the direction of propagation of 
the pump, z, and the components q s ^ (perpendicular to 
the z axis) of their k s ^ vectors are at angles (p s and tpi to 
the x axis. This means that cp is an azimuthal coordinate 
and this will play a central role in the description of the 
orbital angular momentum. The q components can be 
decomposed, on the planes perpendicular to the z axis, 
as 



p s ^i cos (p Syi 
p^sin <p 8ii 




(1) 



where the radial variable, p, extends outwards from z. 
The magnitude of the wave-vector inside the medium 
is ujn/c and for type-I phase matching (eoo) the extra- 
ordinary and ordinary refractive indices are n e (uo p ) and 
no(w s ,i)- Typical values, for example for /3-barium bo- 
rate (BBO), are: at A = 1064nm, n Q = 1.65, n e = 1.54 
and at A = 532nm, n Q = 1.67, n e = 1.55 (26| . 



pump 













idler 

Figure 1: (Color online) Sketch of the SPCD process. A 
gaussian pump propagating in the z direction is incident 

on a short nonlinear crystal. Signal and idler photons 
are produced at angles $ s and di to the pump direction. 

At the output of the nonlinear crystal the two-photon 
state in the wave- vector domain is given by 23|, [27j 



vacuum state and a|(k s ), aj(k^) are creation operators 
for the signal and idler modes with wave vectors k s ,k^, 
respectively. 

Photon pairs generated by parametric down- conversion 
are entangled in arbitrary superpositions of OAM modes 
and we aim to complete this description by including 
the radial behaviour. A natural way to do this is to 
describe the down-converted photons in terms of the 
Laguerre- Gaussian (LG) modes, LG^. Here i corre- 
sponds to the angular momentum carried by the mode, 
£h, and describes the helical structure of the wave front 
around a wave front singularity and p is the number of 
radial zero crossings. 



Calculation of the coincidence amplitudes 

The coincidence probability for finding one signal pho- 
ton in a given LG mode characterized by the i s and p s 
numbers and one idler photon in a given LG mode char- 
acterized by the ti and pi numbers is Pp s s ]p\ = \Cp s s fp\\ •> 
where the coincidence amplitudes Cp*'^. are given by the 
overlap integral 

Cpllpi = (^i^s\lpSPDc) 

= J jdPksdPh^s^) [LG^(k a )]*[LG^(ki)]*(3) 

The pump profile and phase matching condition are given 
by [23 

$(k s ,ki) = [d 3 k p Ep(k p )£(kp - k s - ki)S(uj - w s - Wi) 

(4) 

where the 5 term enforces energy conservation and the 
£ term arises due to phase matching. We can simplify 
this if, instead of the wave vector k, we use its trans- 
verse component q and the corresponding frequency uj. 
As we have assumed a monochromatic gaussian pump of 
frequency uj p = Ui + uo s we can write 



Ep(k p ) = Ep(q)S(uj - uj p ) 

where E p (q) is the Fourier transform of 

[2 1 ( -t 2 
E p(X,(t)) = \ exp 

V TTWp \W£ 



(5) 



(6) 
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the spatial distribution of the pump at the input face of 
the crystal. 

For a crystal of finite thickness, L, in the longitudi- 
nal direction and transverse length much larger than the 
pump beam size, the phase matching condition is [27j 



where $(k s , k^) describes the mode function of the pump 
and the phase matching conditions, |0) is the multimode 
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where Ak z 
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and k z 



^/n 2 (uj)uj 2 I c 2 — q 2 is the longitudinal component of the 
wave vector, k, with transverse component q, angular 
frequency lj and refractive index n(uu). 

If the angle between signal and idler beams is small 
enough that the z-component of the momentum vector 
(\/k 2 — q 2 ) can be approximated by k — q 2 /2k, where 
q = |q|, we can write the phase matching function as 



$(k a ,ki) = E(q s + qi 
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Explicitly: 
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Phase Matching 



It is convenient to calculate the coincidence amplitudes 
in a cylindrical coordinate system, so we re-express (|3]) 
as: 

/ *(p i ,p a ,<p i ,<p a )[LG t p - t {p a ,<p a )Y 
„ Jo 

x [LGp\(pi,<pi)]* pip a dpidp a dtpid<p a , (10) 

where p and ip are the modulus and azimuthal angle, 
respectively, of the transverse component q of the wave 
vector and the normalized LG modes in /c-space are given 

by 
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Here u> is the beam waist (we have assumed z = 0) and 
Lp'(-) is an associated Laguerre polynomial. 

The pump and phase matching functions in (j9]) can 
be written in cylindrical coordinates by performing the 
substitution 



|q* ± Qs| 2 = Pi + p 2 s ± 2^p s cos(^ - <p s ) 



(12) 



to obtain 
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It is then straightforward to see that <&(pi, p s ,(pi,cp s ) 
depends on the radial coordinates in the momentum 
space, pi and p s , and on the difference between the az- 
imuthal angles, (fi — (p 3 . This allows the function to 
be written as a superposition of plane waves with phase 
exp(i£(<pi - cp 3 )): 



$(pi,p s ^i-V s )= J2 h{PuPs)e u ^-*»\ 

£=-oo 

Using this in (JTUJ) , the angular integral becomes 

r 2ir r 2ir 
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which clearly enforces the angular momentum conserva- 
tion, £i + £ a = 0. 

We also re- write the sine function as the inverse Fourier 



Pi + P s ~ 2 PiPs cos((pi -(f s ) \ ^-iL^ 
4k n 
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In this way we calculate 
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(15) where the combinatorial coefficient Kp} }Ps is given by: 
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I, S and T are functions of the dummy variable t: 
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and iF\ is the Gauss hypergeometric function. Although 
this integral is too complicated to be calculated analyti- 
cally, it is simple enough to be evaluated numerically. 

In the limit of a thin crystal, however, we can solve 
the integral analytically (because the integration limits 
depend on the crystal length). This gives 

(i - if + iff' (i + if - if) p ° (-2-rn,) w 
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(20) 



where ji and 7 S are the ratios w p /wi and w p /w Sl so that 
the two 7 factors are the inverse signal and idler widths 
normalized to the pump width. This means that every 
amplitude is invariant under the joint scaling of signal, 
idler and pump widths. 

Comparing (|2Qj) with the numerical evaluation of ([I7|) 
we find excellent agreement for crystal lengths up to tens 
of centimeters, which is two orders of magnitude larger 
than the crystals typically used in experiments. Having 
a thick crystal, therefore, will not greatly affect the mode 
content of the SPDC state, as long as the conditions of 
small emission angles, and similar refractive indices for 
pump, signal and idler fields, which allow Ak z ~ |q^ — 
q s | 2 /2fc p , are met. 



Relation between collinear geometry and thin 
crystals 



where we use the Laguerre-Gaussian modes in the 
real-space form: 
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In this case the integral over the azimuthal coordinate is 
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from which we obtain the well-known conservation law 
for the orbital angular momentum: i s = — ii [291 ]. 

Substituting these into ([22]) we find exactly the same 
expression as ([20]), which means that the amplitudes for 
a collinear beam and a crystal of any experimentally 
realisable size are the same for a quasi-collinear beam in 
the thin crystal approximation. 



STRUCTURE OF THE SPDC STATE 

The advantage of an analytical result over an integral 
one, such as ([T7]h is that it allows us to see more easily 
the role that each parameter plays in determining the 
state of the down-converted photon and the resulting I 
distribution, or spiral bandwidth [23]. Full knowledge 
of this allows the quality of the entangled state to be 
determined and also allows us to compare the width of 
the distribution of t modes of different states. From ([20]) 
it is clear that the distribution of modes is determined 
by the ratio of the signal and idler widths to the pump 
width. We can see this more clearly, and also investigate 
more the effect of taking the radial modes into account, 
if we consider some special cases. 



Note that in collinear geometries (q s = = 0) the 
sine term in the phase matching condition is unity. This 
creates a similarity between having a thin crystal or tun- 
ing a thicker crystal to reach collinear conditions. The 
biphoton state in this case is given by [13, HH 



= J dr$(r)oj(r)o|(r)|0) 



(21) 



where $(r) is the spatial distribution of the pump beam 
at the input face of the crystal, and r is the radial co- 
ordinate in real space. The coincidence amplitudes can 
then be calculated readily from 

Ci s, ti oc / dej) l rdr 



x LG° (r^) [LG l ; s (r^)Y [LG£ 4 M)]* (22) 



Comparison with published results 

Note that for the specific case of p s = pi = eq.([20l 
reduces to 



2\e\ 



(25) 



If we further choose the signal and idler widths to be 
equal, so that 7i = 7 S = 7, then this simplifies to 

(26) 
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This agrees with the result found in [23( as depicted in 
FigH 
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Figure 2: (Color online) Normalized spiral bandwidth 
for 7 S = = 1 and p s = Pi = 0. Results from eq. (10) 
in [23| are in red. 



Although the form of eq. ([26|) is very simple, it is re- 
markably precise when it is compared to the numerical 
evaluation of eq. ([T7|) , where the signal and idler sizes are 
equal and if the eigenstates of pi = p s = have been 
selected. 



Correlation between pi and p s 

It is possible to have an extended view of the structure 
of the entangled system with an array plot that shows 



the contribution PJr 



= \Ci 



for each pair of modes. 



(i) Effect of the pump width on the correlation between pi 
and p s when £ — and signal and idler have the same size. 

We show in Fig|3] the contributions of modes in the 
form \pi,0) (g) \p S: 0). Since 7^ = 7 S in this case, we will 
omit the subscript and simply write 7. On the axes we 
are scanning the discrete values of pi and p s . 

In the limit of an infinite pump width pi and p s are 
delta correlated. Notice how the correlation between pi 
and p s breaks down as the pump beam size approaches 
that of signal and idler. (This corresponds to 7 — » 1). 

In an experiment it is not possible to use a pump 
with arbitrarily large beam waist because the crystal size 
places an effective upper limit. It is also not possible to 
shrink the size of signal and idler enough so that 7 >> 1. 
Therefore, the experimental scenario will be closer to the 
bottom right graph rather than the top left. This means 
that cross correlations between eigenstates of different p 
are to be expected. 




u< 



Figure 3: (Color online) P^° Pi for < pi < 15 and 
< p s < 15, for different pump sizes. The color bar 
shows the color normalisation, which takes place 
between the minimum (0) and the maximum values in 
each graph. The same happens in FigH] and FigJTl 



(ii) Effect of £ on the correlation between pi andp s . 

In the previous section we described the p correlations 
between states of OAM with £ = 0. We now consider the 
case £ 7^ 0. 

It is worth recalling that (in general) the probability 
of selecting a joint state of a given OAM, \£\, decreases 
as the value of \£\ increases. For this reason, if we select 
a specific value of £ : as is done with SLMs, the average 
rate of coincidence counts will be lower than the rate 
measured after selecting states with a lower value of \£\. 
If the value of the selected \£\ is too large, the rate could 
drop too much to obtain a meaningful signal to noise 
ratio and the effects could be masked by noise. This is 
why we show the effect only up to the value \£\ =6. 

A physical reason why the maximum probability of 
detecting modes of similar p shifts to a higher value of p 
for states of a higher OAM eigenvalue can be found by 
considering what happens to the product of the pump 
mode with specific signal and idler modes. In this case 
the key idea is that when we calculate the overlap integral 
between the pump and modes with a given there is 
an optimal combination of p s and pi that maximizes the 
overlap. The higher the value of i, the higher the optimal 
value of the combination of radial indices. 

To show how this happens, consider the product of the 
fundamental pump mode with two pi = p s modes (i.e. 
the integrand of ([TO]) in the short crystal limit). In Fig|5] 
we show three such products for pi = p s = 0, 5, 15, as 
functions of p, omitting the angular dependence. If we 
calculate the (unnormalised) overlap integral, it is pos- 
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Figure 4: (Color online) Correlations between 
< pi < 15 and < p s < 15, for different \£\ values and 
a pump of the same size as signal and idler. Notice how 

the combination of pi and p s with the highest 
probability of being measured shifts along the diagonal 
as higher and higher values of \£\ are selected. Each 
graph is normalised to the maximum value, a 
comparison between the main diagonals for 
\£\ = 4, 5, 6, 7 is given in FigGO 

sible to see (Figj5j) that, although the effect of increas- 
ing the radial indices moves the inner rings towards the 
origin of the coordinates (where the maximum value of 
the gaussian mode of the pump is), the radius of the 
main ring reduces and consequently the overlap will re- 
duce. Thus, as the number of rings increases, an optimal 
value of the overlap integral is reached and subsequently 
it drops, which is exactly what happens along the diago- 
nals of the graphs in FigHJ which are also shown in Fig|6l 



(Hi) Effect of signal-idler size mismatch on the correlation 
between pi and p s . 

In the analysis of the correlations between states of 
different p index we assumed the signal and idler fields 
to have the same size. We now generalize to fields of 
different size. 

When the sizes of signal and idler beams begin to dif- 
fer, the main correlation line between modes of pi and 
p s shifts towards one or the other axis, depending which 
of the beams is larger, thus yielding the highest corre- 
lations between modes of different radial indices rather 
than between modes of equal radial indices. 

FigJTlis a group of four graphs, each featuring the prob- 
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Figure 5: (Color online) Three products in the form 
LG^LGpLGp 1 with £ = 6 and p = 0, 5, 15 (respectively 
in blue, red and brown), plotted radially from the origin 
of the cylindrical coordinate system. 
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Figure 6: (Color online) The overlap between LGq and 

LGpLG~ £ is best achieved at an optimal value of p 
(indicated by red disks), once \£\ is fixed. The red line 
represents the main diagonal in the plot in Fig. 4c and 
the blue line represents the main diagonal in the plot in 
Fig.4d. 



abilities of detecting an eigenstate in the form |p^,0) (g) 
\p s , 0), where the size of the idler beam increases to twice 
the size of the signal beam, which remains of the same 
size of the pump. 

It is clear from the graphs that we can increase the 
detection probability of modes of very different values of 
p simply by changing the relative size of signal and idler 
modes. The consequence of selecting a different £ mode 
is analogous to what is shown in FigHJ 

As we are able to calculate the amplitudes for any val- 
ues of the beams sizes, it is interesting to see what the 
effects of the beams sizes are on the spiral bandwidth. 
We will discuss this in the next section and give multiple 
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(i) Equal radial indices - Equal signal-idler sizes. 

This case is of particular importance because now the 
ratio between pump and signal-idler beams size is the 
only feature that influences the SPDC state. It is a quan- 
tity that is relatively easy to manipulate, for example by 
choosing an initial pump beam size and then magnify- 
ing or de-magnifiying the signal and idler beams. We 
acknowledge that it can be hard to achieve a particular 
magnification, because of limiting apertures in the setup, 
especially between the crystal and the detection instru- 
mentation. Such finite apertures will give a limit on the 
maximum size of the beams, and also losses if this max- 
imum size is exceeded. 

The graphs on the left-hand side of Fig |8] show the SB 
(horizontally) for different values of the ratio of the pump 
width to the signal-idler width, i.e. 7 (which changes on 
the vertical axis). 



Figure 7: (Color online) These graphs show that the 
effect of a width mismatch yields a higher probability of 
finding a state with different radial indices. If, instead, 
the value of 7 S were larger than 7^ the graphs would be 
mirrored with respect to the leading diagonal. 



examples of the effect of varying the beams sizes. 



Effect of beams sizes on the spiral bandwidth. 

We now want to see the effects of the beams sizes on 
the spiral bandwidth (SB). The SB is the collection of 
detection probabilities in the form Pfv^, where the ra- 
dial indices pi and p s have been fixed and we scan the 
£ eigenvalues in an interval (in our case of ±20) around 
£ — 0. Unlike the previous graphs, each graph now does 
not represent the structure of the SPDC state, but rather 
a collection of "slices" of the total modal content, each 
slice representing the SB, i.e. each graph features only 
the £ eigenvalues. 

Each group of graphs will highlight a particular feature 
of the SB determined by the radial indices. Note that the 
normalization has to be performed on each individual 
horizontal line and not on the graph as a whole, because 
the beam size is just a parameter of the detection basis. 
We stress that each value in the graphs corresponds to 
the integral over the whole area where the pump overlaps 
with the detection mode. 

We will divide the effects into four families, each one 
corresponding to equal or different radial indices and 
equal or different signal-idler sizes. In this way we can 
cover all the experimentally interesting effects on the SB. 
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Figure 8: (Color online) These graphs show how the SB 
varies with the pump width, which increases from top 
to bottom of the graphs and its ratio with the signal 
and idler beams is shown on the left. We can see how 
the radial numbers influence this dependence by 
comparing the top with the bottom graphs. The 
positions where the bandwidths that are shown on the 

right are taken are marked on the graphs on the left 
(the blue bandwidth on the right is indicated in yellow 
on the left). 

A larger pump width (or smaller signal and idler 
widths) increases the width of the SB. Also, measuring 
on a basis with larger radial indices gives a larger SB. In 
fact, modes with more rings still give a significant overlap 
with the gaussian pump for larger values of £ than modes 
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with less rings. 

(ii) Equal radial indices - Different signal-idler sizes. 

A rather different and interesting effect is achieved 
when the size of signal and idler differ by some amount. 
Fig|9] shows similar graphs to FigJHJ but here the signal 
and idler beams now have different widths. 




Figure 9: (Color online) These graphs show how the 
difference in the width of signal and idler influences the 
SB. We show the effect for various values of pi = p s . 

A physical explanation of the multiple-branched SBs in 
Fig|9]is given by the fact that the joint detection mode is 
made out of two modes with different widths (therefore 
it will consist of concentric rings of positive and negative 
value) and by the consequence of increasing £: the rings in 
the modes move further away from where the pump mode 
is concentrated, but since they have a different width 
their product will change shape, becoming alternately 
mostly negative or mostly positive. The higher the width 
difference, the higher number of times this process takes 
place in the same amount of £ values. 



The consequence is that the value of the overlap with 
the pump mode (which is always positive) increases and 
decreases alternately as higher and higher t modes are 
chosen, because the mostly positive or mostly negative 
part of the joint detection mode, since they are moving 
away from the center of the beam, will cease to have a 
substantial overlap. 

The number of branches in the graphs in Fig]9] de- 
pends on the relative dimensions of signal and idler, and 
also on the number of rings in the joint detection mode, 
as mostly positive and negative parts can co-exist in the 
same detection mode over many rings. If the pump is too 
small (top of the graphs on the left-hand side, and red SB 
in the graphs on the right-hand side) the effect can't be 
noticed because the gaussian pump always stays mostly 
inside the first ring (so it doesn't matter now many there 
are outside the first one, or their values). However, for a 
larger value of the pump (bottom of the graphs), it ini- 
tially overlaps with many rings, whose overlap can then 
be lost, in the way explained, while we consider higher 
and higher values of The blue SB in the right-hand 
side graphs portrays this effect. 



(Hi) Different radial indices - Equal signal-idler sizes. 

We now consider the case of different values of the 
radial indices pi and p s . Since the probability of detect- 
ing modes of different radial index depends on the ratios 
of the beam sizes, there are choices of widths that may 
significantly increase the probability of detecting some 
particular output modes, as seen in Figd 

In FigJTU] we show the effect of increasing the pump 
width while maintaining the signal and idler at the same 
size, but with different radial indices. 

As long as signal and idler beams retain the same 
width, the number of branches will be always two, in- 
dependent of the values of pi and p s . This behavior is 
due to the fact that the beams have the same width, 
therefore, as higher and higher values of £ are reached, 
the rings in the joint detection mode will be moved away 
from the pump "rigidly" , i.e. maintaining the same shape 
and thus delivering a smooth decay of the overlap. 



(iv) Different radial indices - Different signal-idler sizes. 

We showed that a multiple-branched SB can be ob- 
tained when signal and idler beams have different sizes, 
now we show the same effect while also varying pi and p s . 
To give a more complete description of the possible SBs 
we will now consider the difference between 7^ and j s , 
rather than the ratio. In the first four graphs the ratios 
of the pump width with signal and idler widths have a 
difference (i.e. 7^ — j s ) of 0.5, 0.75, 1 and 2. The pump 
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Figure 10: (Color online) These graphs show how 
different radial indices in the signal and idler modes 
(with equal widths) influence the SB. We show the 
effect for various values of pi , p s . 



width indicated on the left axis of the graphs is relative 
to the signal width. 

The physical explanation for the shapes of the SBs is 
similar to the one supplied in the case of equal radial 
indices and different signal-idler sizes. The effect is de- 
pendent on both the relative size of signal and idler and 
on the number of rings in the modes. The fact that the 
number of branches depends on the smaller of pi and 
p s clarifies what is stated above, namely that it all de- 
pends on how the rings in the modes overlap: if the joint 
detection mode consists of two modes with p s + 1 and 
Pi + 1 rings, the rings will overlap more or less effectively 
depending on both the relative size and the intensity of 
the modes. Therefore, it's the mode with less rings that 
counts. When £ is scanned, the two modes, having a dif- 
ferent width, will shift the rings away from the center 
at different rates, giving alternating mostly positive or 
mostly negative groups of rings. Since the total num- 
ber of rings depends on the lowest of the radial indices, 
the mode with less rings will determine the maximum 
number of branches in the SB, as is shown in FigJTTJ 



CONCLUSION 

We have calculated the form of the SPDC state pro- 
duced with a Gaussian pump beam and investigated the 
resulting correlations in the OAM and radial momentum 
of the down- converted Laguerre- Gaussian beams. Our 
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Figure 11: (Color online) Joint effect of having a signal 
and an idler mode with different radial indices and 
different widths. 



results show excellent agreement with previous works on 
SPDC 0. 

We have shown that, with the usual paraxial and 
collinear assumptions, the phase matching term can be 
neglected allowing an analytic expression for the ampli- 
tudes to be derived. In addition, we have explored the 
effect of the radial contributions and of the beams sizes, 
showing the main families of effects on the radial corre- 
lations and on the SB. 
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In future we will investigate further the phase matching 
term so that we can consider non-collinear systems. We 
also will look at the effect of a spatial light modulator on 
the radial components, since it is the device that is used 
to convert the LG modes for detection by a single mode 
fibre. 
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